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Abstract 

In most of the mass range encompassed by the limits from the direct search and 
the electroweak precision tests, the Higgs boson of the standard model preferably 
decays to bottom quarks. We present, in analytic form, the dominant two-loop 
electroweak correction, of 0(G 2 F m\ ), to the partial width of this decay. It amplifies 
the familiar enhancement due to the OiGprn^) one-loop correction by about +16% 
and thus more than compensates the screening by about —8% through strong- 
interaction effects of order 0(a s GFmf). 
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The standard model (SM) of elementary-particle physics, whose fermion and gauge 
sectors have been impressively confirmed by an enormous wealth of experimental data, 
predicts the existence of a last undiscovered fundamental particle, the Higgs boson H, 
whose mass Mh is a free parameter of the theory. The direct search for the Higgs boson at 
the CERN Large Electron- Positron Collider LEP 2 led to a lower bound of M H > 114 GeV 
at 95% confidence level (CL) [JJ. On the other hand, high-precision measurements, espe- 
cially at LEP and the SLAC Linear Collider SLC, were sensitive to the Higgs-boson mass 
via electroweak radiative corrections, yielding to the value Mh = (85i||^ GeV together 
with an upper limit of Mh < 166 GeV at 95% CL [2J. The vacuum-stability and triviality 
bounds suggest that 130^M^^180 GeV if the SM is valid up to the grand-unification 
scale (for a review, see Ref. [3]). If the Higgs mechanism of spontaneous symmetry break- 
ing, as implemented in the SM, is realized in nature, then we are now being on the eve 
of a groundbreaking discovery, to be made at the CERN Large Hadron Collider (LHC), 
which will go into operation in a just few months from now. After finding a new scalar 
particle, the burning question will be whether it is in fact the Higgs boson of the SM, 
or lives in some extended Higgs sector. Therefore, it is indispensable to know the SM 
predictions for the production and decay rates of the SM Higgs boson with high preci- 
sion. Its decay to a bottom-quark pair, H — > bb, is of paramount interest, as it is by far 
the dominant decay channel for M^^140 GeV (see, e.g., Ref. [1]). On the other hand, 
the inverse process, bb — > H, was identified to be a crucial hadroproduction mechanism, 
appreciably enhancing the yield due to gluon fusion [5]. Precise knowledge of the bottom 
Yukawa coupling is also requisite for reliable predictions of associated hadroproduction of 
Higgs bosons and bottom quarks [6J. 

The purpose of this Letter is to fill a long-standing gap in our knowledge of the 
quantum corrections to the partial width r& of the H — > bb decay, by providing, in analytic 
form, the dominant two-loop electroweak correction, of 0(G 2 F mf), where Gf is Fermi's 
constant and mt is the top-quark mass. This correction also applies to the cross section 
of bb — > H. Surprisingly, it turns out to be more than twice as large as the 0{a s GF r m 2 ) 
one, which is formally enhanced by one power of the strong-coupling constant a s . In 
the discussion of virtual top-quark effects, it is useful to distinguish between universal 
corrections, which are independent of the produced fermion flavor, and non-universal 
corrections, which are specific for the H — > bb decay because bottom is the weak-isospin 
partner of top. Here, we have to consider both types. 

Prior to going into details with our calculation, we briefly review the current status of 
the radiative corrections to Y b in the intermediate mass range, defined by M w < M H < 
2Mw- As for effects arising solely from quantum chromodynamics (QCD), the full m& 
dependence is known in 0(a s ) [7\. In 0(a 2 ), the leading [8] and next-to-leading [9] terms 
of the expansion in m 2 /M H of the Feynman diagrams without top quarks are available. 
Those involving top quarks either contain gluon self-energy insertions or represent cuts 
through three-loop double-triangle diagrams; the former contribution is exactly known 
[TO] , while the four leading terms of the expansion in are known in the latter 

case [llj. In O{o i s ), the diagrams containing only light degrees of freedom were evaluated 
directly [12J, while those involving the top quark were treated in the framework of an 
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appropriate effective field theory [13] . As for purely electroweak corrections, the one- 
loop result is completely known [T3]. At two loops, the dominant universal correction, 
of 0(Gj?mf), was already studied in Ref. [15J, while the non-universal one is considered 
here for the first time. As for mixed corrections, the universal [16] and non-universal p2] 
0(a s GFm^) terms at two loops and the universal [32] and non-universal [19J O^G^m^) 
terms at three loops are available. 

We now outline the course of our calculation and exhibit the structure of our results. 
Full details will be presented in a forthcoming communication [20]. For convenience, 
we work in 't Hooft-Feynman gauge. As usual, we extract the ultraviolet divergences 
by means of dimensional regularization, with D = 4 — 2e space-time dimensions and 
't Hooft mass scale /i. We do not encounter ambiguities related to the treatment of 75 
in D dimensions and are thus entitled to use the anti-commuting definition. We adopt 
Sirlin's formulation of the electroweak on-shell renormalization scheme [21], which uses Gf 
and the physical particle masses as basic parameters. We take the Cabibbo-Kobayashi- 
Maskawa quark mixing matrix to be unity, which is well justified because the third quark 
generation is, to good approximation, decoupled from the first two [22] . For convenience, 
we renormalize the Higgs sector by introducing counterterm vertices involving tadpole and 
Higgs-boson mass counterterms, 5t and 5Mh, respectively [23]. Specifically, 5t is adjusted 
so that it exactly cancels the sum of all one-particle-irreducible tadpole diagrams. 

Detailed inspection reveals that, to the orders considered here, the amputated matrix 
element of H — ► bb exhibits the simple structure 

A = A + B{f-f)u_, (1) 

where u± = (1 ±7s)/2 are the helicity projection operators, p and p are the four-momenta 
of b and b, respectively, and A and B are Lorentz scalars. Including the wave-function 
renormalizations of the external particles and employing the Dirac equation, we find the 
transition matrix element to be 

T = {z~ H (y]Z hiL Z h , R A + m b Z h . L B) s, (2) 

where s = u(p,r)v(p,r), with r and r being spin labels. Owing to parity violation, the 
left- and right-handed components of the bottom-quark field, b L R = u T b, participate dif- 
ferently in the electroweak interactions and thus receive different wave-function renormal- 
izations, Zb t L/R- At tree-level, we have A^ = —m^jv and = 0, where v = 2~ 1 / 4 G F 1 ^ 2 
is the Higgs vacuum expectation value. Here and in the following, superscripts enclosed 
in parentheses denote the loop order. In Sirlin's formulation of the electroweak on-shell 
scheme, where Fermi's constant is introduced to the SM through a charged-current pro- 
cess, namely muon decay, the SU(2) gauge coupling g = 2Mw/v does not receive power 
corrections in mt, so that [21] 

M Wfi _ M w 



to the orders considered here, which implies that the renormalization of v is reduced to 
the one of Mw- Here and in the following, bare quantities carry the subscript 0. It hence 
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follows that we need to perform a genuine two-loop renormalization of Zh, m b , Z b ,L/R, and 
Myy, while a one-loop renormalization of M# and m t is sufficient. As usual, we denote 
the sums of all one-particle-irreducible H, f (f = b,t), and W self-energy diagrams at 
four-momentum transfer q as z£#(g 2 ), i[d (g 2 ) + u+Ef^q 2 )) + m/^E/.st^ 2 )], and 

-i\(g^ v - q^q v /q 2 )Y, W)T (q 2 ) + [<fq v j 'q 2 )Z W)L (q 2 )}, and split the bare masses as M 2 H/W0 = 
M]j! W + 5Mjj/ w and m/ j0 = rrif + 8m f. Imposing the on-shell renormalization conditions 
on the dressed propagators then yields 
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where 



(i + £/,L(g 2 ))(i + £/,*(g 2 )) 

Relations that, to the order of our analysis, are equivalent to Eq. (jlj) were found in Ref. 
using an alternative approach. 

Performing a loop expansion and eliminating all bare masses, we thus obtain 
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The Feynman diagrams contributing to and Bq Z> are depicted in Fig. [TJ They 
are generated and drawn using the program FeynArts [26] and evaluated using the pro- 
gram MATAD [27J, which is written in the programming language FORM |28j, by applying 
the asymptotic-expansion technique (for a careful introduction, see Ref. [29]). Here, \ 
and <p denote the neutral and charged Higgs-Kibble ghosts with masses Mz and Mw, 
respectively. The crosses in Figs. [TJs) and (t) indicate the insertions of the Higgs-boson 
mass and tadpole counterterms iSt/vo and —i (8t/vo + 5Mfj) in a (/>-boson line and 
a Hcjxp vertex, respectively. In the soft-Higgs limit, M# <C mt, which is underlying our 
analysis, the diagrams in Figs. [T^a)-(s) can also be evaluated by applying a low-energy 
theorem (see Ref. [SU] and references cited therein) to the corresponding 6-quark self- 
energy diagrams that emerge by removing the external Higgs-boson line. This provides 
a powerful check for our calculation. Apart from the diagrams in Fig. [TJ we also need to 
calculate the relevant one-particle-irreducible H, b, and W self-energy diagrams at two 
loops. Furthermore, we need to expand all the relevant one-loop diagrams through 0(e). 
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We are now in a position to present our final results for the universal correction 
parameter 5 U and the relative correction to r&. They read 

7 
6 
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where N c = 3 and C F = (iV 6 2 - 1)/(2N C ) = 4/3 are color factors, x t = (G F m 2 t ) / '{8tt 2 \/2) , 
C(2) = tt 2 /6, and 

r ( 0) = V2N c G F M H mj ^ _ 4m|y /2 _ 



8^ 



If we convert Eq. ([9]) to a mixed renormalization scheme which uses the on-shell definitions 
for the particle masses and the definitions of the modified minimal-subtraction (MS) 
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Figure 1: Diagrams contributing to H — > bb at 0(G 2 F mf). 

scheme for all other basic parameters, then we find agreement with Eq. (15) for x = in 
the paper by Djouadi et al. |15j . However, the corresponding result for the electroweak 
on-shell scheme presented in their Eq. (27) for x = disagrees with our Eq. (jUJ). We 
can trace this discrepancy to the absence in their Eq. (25) of the additional finite term 
(j^ApW which arises from the renormalization of the one-loop result in their Eq. (7) 
according to the prescription in their Eq. (18). The 0(G 2 F mf) term in Eq. ( ITUl) represents 
a new result. 

In Eqs. dU) and tTTU]) . we have also included the two-loop C(a s Gi?mf) corrections 
[T6][T7] . which we reproduced using our calculational techniques. As for the QCD renor- 
malization, it is understood that appearing in Eq. (TTTj) is defined in the MS scheme 
as m;, = m&(Mff), while the electroweak part of the renormalization remains in the on- 
shell scheme. This modification ensures that large logarithms of the type In (Mfj/m 2 ,) 
that would otherwise appear already at 0(a s ) and spoil the convergence behavior of the 
perturbation expansion are properly resummed according to the renormalization group 
(RG) [7J. Since we wish to treat m t on the same footing as m&, we adopt this mixed 
scheme for m t as well. The analysis at O^alGprn 2 ) [18119] reveals that Eqs. (jHJ) and ffTOl 
may be further RG-improved by taking m t and a s to be m t = m t (m t ) and a s = af\mt), 
respectively. 

Finally, we explore the phenomenological implications of our results. Adopting from 
Ref. [22] the values G F = 1.16637 x 10" 5 GeV" 2 , af\M z ) = 0.1176, M z = 91.1876 GeV, 
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Table 1: Relative corrections to T T , T c , and Tb at 0(G F m 2 ), 0(G%mf), and 0(a s G F m 2 ). 



Order 


rv/if) r c /r(°) ryrf 


0{G F m 2 ) 
0{G\m\) 
0(a s G F m 2 ) 


+2.021% +2.021% +0.289% 
+0.064% +0.064% +0.047% 
+0.060% +0.452% -0.022% 



and mf ole = 174.2 GeV for our input parameters, so that af\m t ) = 0.1076 and m t = 
166.2 GeV, we evaluate Eqs. © and ([TOD to 0(G F m 2 ), 0(G 2 F mj), and 0(a s G F m 2 ). For 
comparison, we also evaluate the relative corrections to Ti and T q , where I = e,fi,T and 
q = u,d, s, c, which, to the orders considered here, are given by 

-§) = (i+<y 2 -i, (12) 
1 1 

-^ = (l + A QC D)(l + <y 2 -l, (13) 
1 g 

where [7] 

A Q cd = -rC F - (14) 

7T 4 

is the 0(a s ) correction in the limit m q <C M#. 

The results are listed in Table [U We observe that the 0(G F mf) correction to 
increases the enhancement due to the 0(G F m 2 ) one by about 16% and has more than 
twice the magnitude of the negative 0(a s G F m 2 ) one. Also in the case of Ti, the 0(Gpmj) 
correction exceeds the 0(a s G F m 2 ) one. The situation is quite different for the case of 
T q , which is due to the additional appearance of the sizeable product term 2Aq CD 5^ 1 ^ in 
Eq. pp. 

In conclusion, we analytically calculated the dominant electroweak two-loop correc- 
tion, of order 0(G 2 F m^), to the H — > bb decay width T b of an intermediate-mass Higgs 
boson, with Mh *C m t . We performed various checks for our analysis. The ultraviolet di- 
vergences cancelled through genuine two-loop renormalization. Our final result is devoid 
of infrared divergences related to infinitesimal scalar-boson masses. We reproduced those 
Hbb vertex diagrams where the external Higgs boson is coupled to an internal top-quark 
line, which we had computed directly, through application of a low-energy theorem. After 
switching to a hybrid renormalization scheme, our 0(G F mf) result for the universal cor- 
rection 5 U agrees with Ref. [TB] . Using our computational techniques, we also recovered 
the 0(a s G F m 2 ) corrections to 5 U and T b . The 0(G 2 F mf) correction to T b amplifies the 
familiar enhancement due to the 0(G F m 2 ) correction by about +16% and thus more than 
compensates the screening by about —8% through QCD effects of 0(a s G F m 2 ). 
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